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We present a measure of entanglement that can be com- 
puted effectively for any mixed state of an arbitrary bipartite 
system. We show that it does not increase under local manip- 
ulations of the system, and use it to obtain a bound on the 
teleportation capacity and on the distillable entanglement of 
mixed states. 



I. INTRODUCTION 

In recent years it has been realized that quantum 
mechanics offers unexpected possibilities in information 
transmission and processing, and that quantum entangle- 
ment of composite systems plays a major role in many 
of them. Since then, a remarkable theoretical effort has 
been devoted both to classifying and quantifying entan- 
glement. 

Pure-state entanglement of a bipartite system is 
presently well-understood, in that the relevant param- 
eters for its optimal manipulation under local operations 
and classical communication (LOCC) have been identi- 
fied, in some asymptotic sense |l) as well as for the single 
copy case [^. Given an arbitrary bipartite pure state 
\iPab}, the entropy of entanglement E{^ab) namely 
the von-Neumann entropy of the reduced density ma- 
trix Pa = Tt:b\iPab){''Pab\, tells us exhaustively about 
the possibilities of transforming, using LOCC, \iPab) into 
other pure states, in an asymptotic sense. When manip- 
ulating a single copy of this information is pro- 
vided by the n entanglement monotones Ei = X)r=; 
{I = l,...,n) |], where are the eigenvalues of pA in 
decreasing order. 

Many efforts have also been devoted to the study of 
mixed-state entanglement. In this case several measures 
have been proposed. The entanglement of formation 
Ef{p) §1 — or, more precisely, its renormalized version, 
the entanglement cost Ec{p) ||] — and the distillable en- 
tanglement Eu{p) 1^ quantify, respectively, the asymp- 
totic pure-state entanglement required to create p, and 
that which can be extracted from p, by means of LOCC. 
The relative entropy of entanglement ||] appears as a 
third, related measure Q which interpolates between Ec 
and Ed §. 

However, in practice it is not known how to effectively 
compute these measures, nor any other, for a generic 



mixed state, because they involve variational expres- 
sions. To our knowledge, the only exceptions are Woot- 
ter's closed expression for the entanglement of formation 
Ep{p) [and concurrence C {p)] oftwo-qubit states [pf , and 
its single-copy analog E2 (p) also for two qubits |9||^ 

Multipartite pure-state entanglement represents the 
next order of complexity in the study of entanglement, 
and is of interest, because one hopes gain a better under- 
standing of the correlations between different registers of 
a quantum computer. Consider a tripartite state ItJjabc) ■ 
Some of its entanglement properties depend on those of 
the two-party reduced density matrices, which arc in a 
mixed state. For instance, the relative entropy of pab = 
trc\''pABc){''pABc\ has been used to prove that bipartite 
and tripartite pure-state entanglement are asymptoti- 
cally inequivalent Thus, the lack of an entangle- 
ment measure that can be easily computed for bipartite 
mixed states is not only a serious drawback in the study 
of mixed-state entanglement, but also a limitation for 
understanding multipartite pure-state entanglement. 

The aim of this paper is to introduce a computable 
measure of entanglement , and thereby fill an impor- 
tant gap in the study of entanglement. It is based on the 
trace norm of the partial transpose p'^^ of the bipartite 
mixed state p, a quantity whose evaluation is completely 
straightforward using standard linear algebra packages. 
It essentially measures the degree to which p^'^ fails to 
be positive, and therefore it can be regarded as a quan- 
titative version of Peres' criterion for separability [p^ . 
From the trace norm of p^^, denoted by Hp"^"*]]!, we will 
actually construct two useful quantities. The first one is 
the negativity, 



Af{p) 



(1) 



which corresponds to the absolute value of the sum of 
negative eigenvalues of p"^-* and which vanishes for 
unentangled states. As we will prove here, M{p) does not 
increase under LOCC, i.e., it is an entanglement mono- 



tone 1 14 1, and as such it can be used to quantify the 
degree of entanglement in composite systems. We will 
also consider the logarithmic negativity, 



Eu{p) EE iog2 Hp' 



(2) 



which again exhibits some form of monotonicity under 
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LOCC (it does not increase during deterministic distilla- 
tion protocols) and is, remarkably, an additive quantity. 

The importance of M and Ef^ is boosted, however, be- 
yond their practical computability by two results that 
link theses measures with relevant parameters charac- 
terizing entangled mixed states. The negativity will be 
shown to bound the extent to which a single copy of 
the state p can be used, together with LOCC, to perform 
quantum teleportation ||l5|. In turn, the logarithmic neg- 
ativity bounds the distillable entanglement E'jj contained 
in p, that is, the amount of "almost pure" -state entan- 
glement that can be asymptotically distilled from p**^, 
where "almost" means that some small degree e of imper- 
fection is allowed in the output of the distillation process. 

Remarkably, this last result has already found an ap- 
plication in the context of asymptotic transformations 
of bipartite entanglement ||l6|| , as a means to prove that 
(PPT) bound entangled states jl^ cannot be distilled 
into entangled pure states even if loaned (i.e. subse- 
quently recovered for replacement) pure-state entangle- 
ment is used to assist the distillation process. In this way, 
the bound on distillability implied by Ej\f has contributed 
to prove that, in a bipartite setting, asymptotic local ma- 
nipulation of mixed-state entanglement is sometimes, in 
contrast to its pure-state counterpart, an inherently irre- 
versible process. 

We have divided this paper into VII sections. In sec- 
tion II some properties of the negativity A/", such as its 
monotonicity under LOCC, and of the logarithmic neg- 
ativity Ej^ are proved. We also discuss a more general 
construction leading to several other (non-increasing un- 
der LOCC) negativities. In sections III and IV we derive, 
respectively, the bounds on teleportation capacity and on 
asymptotic distillability. Then in section V we calculate 
the explicit expression of M and £V for pure states and 
for some highly symmetric mixed states, also for Gaus- 
sian states of light field. In section VI extensions of these 
quantities to multipartite systems are briefly considered, 
and section VII contains some discussion and conclusions. 

II. MONOTONICITY OF M(p) UNDER LOCC. 

In this section we show that the negativity M{p) is an 
entanglement monotone. We first give a rather detailed 
proof of this result. Then we sketch an argument extend- 
ing this observation to several other similarly constructed 
negativities — like the robustness of entanglement p^ . 

A. Definition and basic properties 

From now on we will denote by p a generic state of 
a bipartite system with finite dimensional Hilbert space 
Ua®'Hb = C-* (g) C''^ shared by two parties, Alice 
and Bob. p-^'* denotes the partial transpose of p with 



respect to Alice's subsystem, that is the hermitian, trace- 
normalized operator defined to have matrix elements, 

= {kA,3B\p\iA,lB), (3) 

for a fixed but otherwise arbitrary orthonormal product 
basis |m, js) = ® \j) ^ G Ha ® Ti-B- The trace 
norm of any hermitian operator A is \\A\\i = trV A'^A 
[]l9| . Sec. VI. 6], which is equal to the sum of the absolute 
values of the eigenvalues of A, when A is hermitian [^ . 
For density matrices all eigenvalues are positive and thus 
IIpIIi = trp = 1. The partial transpose p^^ also satisfies 
tr[p^-*] = 1, but since it may have negative eigenvalues 
Pi < 0, its trace norm reads in general 

||p^-|li = l + 2|^M.|^l + 2AA(p). (4) 

i 

Therefore the negativity Af{p) — the sum | J^i I of the 
negative eigenvalues pi of p^^ — measures by how much 
p^^ fails to be positive definite. Notice that for any sep- 
arable or unentangled state ps PH , 

Ps = ^Pk\ekJk){ekJk\; Pfc>0,^Pfc = l, (5) 

fc fc 

its partial transposition is also a separable state [|2|, 

P^'^E^'fcl^fe'/'^X^'/'^I^O, (6) 
fc 

and therefore H/of-^Hi = 1 and Af{ps) = 0. 

The practical computation of J^{p) is straightforward, 
using standard linear algebra packages for eigenvalue 
computation of hermitian matrices. On the other hand, 
this representation is not necessarily the best for proving 
estimates and general properties of M{p). To begin with 
a simple example, consider the property that N{p) does 
not increase under mixing: 

Proposition 1: M is a convex Junction, i.e., 

n{^p..p^ <Y.P.N{p.) , (7) 

whenever the pi are hermitian, and pi > with pi = 
1. 

There is nothing to prove here, when we write JV{p) = 
(IIp"^"*!!! ~ l)/^! and observe that || • ||i, like any norm, 
satisfies the triangle inequality and is homogeneous of 
degree 1 for positive factors, hence convex. 

However, the fact that ||p||i is indeed a norm is not so 
obvious, when it is defined in terms of the eigenvalues. 
This is shown best by rewriting it as a variational expres- 
sion. Our reason for recalling this standard observation 
from the theory of the trace norm is that the same varia- 
tional expression will be crucial for showing monotonicity 
under LOCC operations. The variational expression is 
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simply the representation of a general hermitian matrix 
A as a difference of positive operators: Since we are in 
finite dimension we can always write 



A 



a+p' 



(8) 



where > are density matrices (tr[p^] = 1) and a± > 
are positive numbers. Note that by taking the trace of 
this equation we simply have tr[A] = a+ — a_. 

Lemma 2: For any hermitian matrix A there is a de- 
composition of the form ^ for which a+ + a_ is minimal. 
For this decomposition, ||A||i = a+ + a_, and a_ is the 
absolute sum of the negative eigenvalues of A. 

Proof: Let be the projector onto the negative eigen- 
valucd subspace of A, and J\f = —tr[AP~] the absolute 
sum of the negative eigenvalues. We can reverse the 
decomposition (|^) to obtain that A + a^p^ is positive 
semidefinite. This implies that 

< tr[(^ + a^p-)p-] = -7V+ a_tr[p-p-] (9) 

But tr[p^P^] < 1, that is a_ > TV. This bound can 
be saturated with the choice a_p~ = —P~AP~ (corre- 
sponding to the Jordan decomposition of A, where p~ 
and p+ have disjoint support), which ends the proof. □ 
For the negativity we therefore get the formula 



Af{A) ^ inf |a_ 



a^p 



(10) 



where the infimum is over all density matrices p^ and 
a± > 0. 

Another remarkable property of J^{p) is the easy way 
in which M{pi ® P2) relates to the negativity of pi and 
that of p2- This relationship is an important, but notori- 
ously difficult issue for discussing asymptotic properties 
of entanglement measures (see, e.g., |^ for a discussion 
and a counterexample to the conjectured additivity of 
the relative entropy of entanglement). 

For the entanglement measure proposed in this paper 
we get additivity for free. We start from the identity 
\\pi ® /O2II1 = IIpiIIi IIP2II1, which is best shown by us- 
ing the definition of the trace norm via eigenvalues, and 
observe that partial transposition commutes with taking 
tensor products. After taking logarithms, we find for the 
logarithmic negativity: 



Em{Pi ® P2) = Ej^{pi) + Em{p2) 



(11) 



It might seem from this that £\r is a candidate for 
the much sought for canonical measure of entanglement. 
However, it has other drawbacks. For instance, it is not 
convex, as is already suggested by the combination of 
a convex functional (the trace norm) with the concave 
log function, which implies that it increases under some 
LOCC. And although it has an interesting, monotonic 
behavior during asymptotic distillation (as shown in sec- 
tion IV) , it does not correspond to the entropy of entan- 
glement for pure states (see section V). 



B. Negativity as a mixed-state entanglement 
monotone. 

By definition, a LOCC operation (possibly for many 
parties) consists of a sequence of steps, in each of which 
one of the parties performs a local measurement, and 
broadcasts the result to all other parties. In each round 
the local measurement chosen is allowed to depend on the 
results of all prior measurements. If at the end of a LOCC 
operation with initial state p the classical information 
available is "i", which occurs with probability pi, and 
final state conditional on this occurrence is p'^, we require 
of an entanglement monotone |14| E that 



e{p)>Y^pM) 



(12) 



It is clear by iteration that this may be proved by look- 
ing at just one round of a LOCC protocol, consisting of 
a single local operation. In the present case, since N 
makes no distinction between Alice and Bob, it suffices 
to consider just one local measurement by Bob. 

Now the most general local measurement is described 
by a family AAi of completely positive linear maps such 
that, in the notation used in the previous paragraph, 
Mi{p) = pip[. These maps satisfy the normalization con- 
dition tr(A^i(p)) = tr(p). This can be further simpli- 
fied when some M.i can be decomposed further into 
completely positive maps, e.g., Mi = M.[ +M'/. Then 
we may simply consider the finer decomposition as a finer 
measurement, with the result i replaced by two others, 
i' and i" . Using the convexity already established it is 
clear that it suffices to prove (|l^ ) for the finer measure- 
ment. That is, we can assume that there are no proper 
decompositions of the A4i, or that A4i is "pure". This 
is equivalent to Mi taking pure states to pure states, or 
to the property |23| that it can be written with a single 
Kraus summand. Taking into account that this describes 
a local measurement by Bob, we can write 



M^{p) = [Ia ® M,) p {Ia ® m}) 



(13) 



where the Kraus operators Mi must satisfy the normal- 
ization condition M^Mi < Ib- For computing the 
right hand side of Eq. dl3) we need that 



M.,{pf^ ^ M^i^P^-') 



(14) 



which immediately follows from (|l^) by expanding p as a 
sum of (not necessarily positive) tensor products. A sim- 
ilar formula holds for Alice's local operations, but with a 
modified operation Mi on the right hand side, in which 
the Kraus operators have been replaced by their complex 
conjugates. 

Consider the decomposition 



pTA = {l+N)p+ - Np- 



(15) 
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with density operators and N = Af{p) . Then we can 
also decompose the partially transposed output states: 

= il + N)M^{p+) ~ NM,{p-). (16) 

Dividing by pi we get a decomposition of precisely the 
sort ( p^ defining J\f{p[). The coefficient a_ — N/pi 
must be larger than the infimum, i.e., M{p[) < N/pi. 
Multiplying by pi and summing, wc find the following 
inequality. 

Proposition 3: 

Y,p,N[p[)<N{p), (17) 

i 

i.e., M{p) is indeed an entanglement monotone. 



C. Other negativities. 

Both the proof of convexity and the proof of mono- 
tonicity are based on the variational representation of 
the trace norm in Lemma 2. The abstract version of 
this Lemma is the definition of the so-called base norm 
II • lis associated with a compact set S* in a real vector 
space The negativity introduced above then corre- 
sponds to a special choice of S", and we can easily find the 
property of S required for proving LOCC monotonicity 
in the abstract setting. Other choices of S then lead to 
other entanglement monotones, some of which have been 
proposed in the literature. 

For our purposes we can take S as an arbitrary com- 
pact convex subset of the hermitian operators with unit 
trace, whose real linear hull equals all hermitian opera- 
tors. Then, in analogy to Lemma 2, we define the asso- 
ciated base norm and "S'-negativity" as 



P||5-inf{c 
Ns{A) = inf {c 



A = a+p+ - a_p", (18) 
a± > 0, e 5*1, 
A = a+p+-a_p", (19) 
a± > 0, p± e S*}, 

(20) 



Note that once again, if A has trace 1 we have that 
ll^lls = 1 + 2Ns{A). Then norm and convexity prop- 
erties of Ms and || • \\s follow exactly as before. 

Taking S as the set of all density matrices, we get 
\\A\\s = ll^lli, for all hermitian A, and a totally unin- 
teresting entanglement quantity, as Ms{p) vanishes for 
all density matrices. The negativity of the previous sec- 
tion corresponds to the choice of S equal to the set of all 
matrices A such that A = A\ tvA = 1, and A^^^ > (ad- 
ditionally, we have replaced A^^ with A in the LHS of eq. 
( p^ ) A, so that we can write Af{p) instead of N{p'^^)). 



We could have also taken S as the subset of density 
matrices with positive partial transpose, /c^ > and 
p^-^-* > 0. In this case S corresponds to all states such 
that its partial transpose is also a state. The resulting 
quantity we will denote by Mppts- Even more restric- 
tively, if we take for S the set of separable density oper- 
ators, i.e. we take p^ (and therefore also p^"^'*) in eqs. 
( p^ , p^ to be separable, the corresponding quantity Nss 
amounts to the robustness of entanglement, originally in- 
troduced in |l^ (see also |^^) as the minimal amount 
of separable noise needed to destroy the entanglement of 
p. From the inclusions between the respective sets S we 
immediately get the inequalities 



(21) 



In general, all these inequalities are strict. For exam- 
ple, Mss{p) vanishes only on separable states, whereas 
J^pptsip) and Af{p) vanish for all ppt-states. 

We claim that also Afss and Afppts are entanglement 
monotones. The proof is quite simple. An analysis of 
the arguments given in the previous section shows that 
we really used only one property of S, namely that for 
all operations Mi appearing in a LOCC protocol, we 
have Mi{p) S , whenever p G , where notes 
the cone generated by S (equivalently the set of Xp with 
X > 0, p G S). But this is obvious for both separable 
states and ppt-states. 



III. UPPER BOUND TO TELEPORTATION 
CAPACITY 

Sections III and IV are devoted to discuss applica- 
tions of the previous results. More specifically, we de- 
rive bounds to some properties characterizing the entan- 
glement both of a single copy of a mixed state p (this 
section) and of asymptotically many copies of it (next 
section). 

For a single copy of a bipartite state p acting on 
C*^ ® C''^, where we set di — d2 = m for simplicity, 
an important question in quantum information theory 
is to what extend this state can be used to implement 
some given tasks requiring entanglement, such as telepor- 
tation. The best approximation Popt{p) to a maximally 
entangled state. 



I* 



1 



\aA <i$ as), 



(22) 



that can be obtained from p by means of LOCC is then 
interesting, because it determines, for instance, how use- 
ful the state p is to approximately teleport log2 m qubits 
of information. In this section we will show that the neg- 
ativity J^if) provides us with an explicit lower bound 
on how close p can be taken, by means of LOCC, to the 
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state From here a lower bound on the teleportation 
distance (i.e., an upper bound on how good teleportation 
results from p) will also follow. 

A. Singlet distance 

In order to characterize the optimal state Popt{p) 
achievable from p by means of LOCC, we need to quan- 
tify its closeness to the maximally entangled state P+ = 
|$+)($+|. Let pi and p2 be two density matrices. The 
trace norm of pi— P2, (or absolute distance |^^), is a mea- 
sure of the degree of distinguishability of pi and p2, and 
it is therefore reasonable to use it to measure how much 
P{p) — the state resulting from applying a local proto- 
col P to state p — resembles P+. In what follows we will 
prove that the negativity is a lower bound to the singlet 
distance of p, 



A(P+,p) 



inf I 

p 



^(P)lli 



(23) 



where the infimum is taken over local protocols P. 

We start be recalling that the absolute distance 
D{pi,p2) = IIpi — P2II1 is a convex function ||2^ , 



^PiD{a,pi) > D{cr,^piPi), 



(24) 



We can then assume that Popt{p) has already under- 
gone a twirling operation. This means that it is a noisy 
singlet pT], 



Pp =pP+ + (l-p)^-^, 



(27) 



from which the absolute distance to P+ can be easily 
computed, D{P+,pp) — 2(1 — p){m? — l)/m?. Similarly, 
the trace norm of p^^ reads — mp + (1 — p)/m, 

and therefore 



DiP+,pp)^ 2(1 



-)■ 



(28) 



The lower bound to the singlet distance ( 23 ) follows now 
straightforwardly from the monotonicity of ll/C"^"*!!! (or 



M{p)) under LOCC, that is, > 
and therefore 



\Poptip) 



A(P+,p)>2(l 



\\P 



m 



(29) 



Therefore we have proved the following bound for the 
singlet distance. 
Proposition 4: 



A(P+,p)>2(l- 



1 + 2J\f{p) 



m 



)• 



(30) 



which confirms, as already assumed, that the optimal ap- 
proximation P{p) to P+ can always be chosen to be a sin- 
gle state — as opposed to a distribution of states {pi,Pi} 
corresponding to the output of a probabilistic transfor- 
mation. Therefore, in eq. ( p3| ) we need only consider 
deterministic protocols P based on LOCC. 

A second feature of the absolute distance that we need 
is that 



D{WpiW\ Wp2W^) = D{pi,p2) 



(25) 



for any unitary transformation W. Properties (|2J) and 
( p5| ) together imply that the best approximation to the 
maximally entangled state P+ can always be "twirled" 
without losing optimality. Consider the state 



dU U ® U* Poptip)U^ ® U'^* , 



(26) 



which the parties can locally obtain from Popt{p) by Al- 
ice applying an arbitrary unitary U , by Bob applying U* ^ 
and by then deleting the classical information concern- 
ing which unitary has been applied. It follows from the 
invariance of P+ under U ®U* and from property ( |25| ) 
that D[U®U*Popt[p)U^®U^\P+) = D{Popt{p).P+) for 
any U . Then property (|2^) implies that the mixture in 
eq. ( |26| ) is not further away from P+ than Popt(p). But 



Ppvtip) was already minimizing eq. (P3[), and therefore 



opt 

state 



must also be optimal. 



B. Teleportation distance 

A quantum state p shared by Alice and Bob can be 
used as a teleportation channel A ||l5|. That is, given 
the shared state p and a classical channel between the 
parties, Alice can transmit an arbitrary (unknown) state 
(f> e C" to Bob with some degree of approximation. Let 
^T,p{4') be the state that Bob obtains when Alice sends 4> 
using p and some protocol T involving LOCC only. The 
teleportation distance 



d{K)= / d^D{cp,K{cp)), 



(31) 



where D{<p,A{<p)) = \\\(t>){(t>\ - H(t>)\\i, can be used to 
quantify the degree of performance of the channel. The 
measure d(j) is consistent with the Haar measure dU in 
SU{m), and thus d{A) is invariant under the twirling 
of the channel, that is the application of an arbitrary 
unitary U to 4> previous to the teleportation, followed 
by the application of [/^ after the teleportation scheme. 
Indeed, 



d(A) = dW D{WPoW\A(WPoW'')), 



(32) 



for some reference state Pq 
(Eq) of the trace norm, eq. 



|0o)((/)o|, and using property 
is also equal to 
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J dW d{WPoW'',U''A{UWPoW''U^)U). (33) 

We can now average over U to obtain 

d{A)^JdU J d(f>D{(f>,U^A{U(t))U), (34) 

where the right side of the equation corresponds to the 
teleportation distance of the twirled channel. 

We next adapt a reasoning of the Horodeckis |2^] to 
our present situation. It uses an isomorphism between 
states PA and channels A due to Jamiolkowski Let 
us ascribe the channel A to the state pA = {I (g) A)P+. 
The state pA can be produced by sending Bob's part 
of the bipartite system in state down the channel 
A. Conversely, the standard teleportation protocol fl^ j 
(or a slight and obvious modification of it) applied to 
state p\ reproduces the channel A with probability l/m? . 
However, if the state p^ is a noisy singlet pp, then the 
corresponding channel is the depolarizing channel 

A'i'^{Q)=PQ+{l-p)^, (35) 

which the standard teleportation scheme reproduces with 
certainty using state pp. For this case d{K'^'^P) = 2(1 — 
p)(m — l)/m. Therefore there is a complete physical 
equivalence between noisy singlets and depolarizing tele- 
portation channels. In addition, 

rf(Af^) = ^i?(P+,Pp). (36) 

Now, since both quantities d and D are invariant under 
twirling, and any channel (state) can be taken into the 
depolarizing (noisy singlet) form, this equality holds for 
any channel A and state pA- 

Lemma 5 (adapted from ||2^): The minimal distance 
dmin{p) that can be achieved when using the bipartite 
state p to construct an arbitrary teleportation channel is 
given by 

Tfi 

dmm(p) = ^-A(P+,p). (37) 

m + 1 

Proof. dmin{p) < mA{P+,p)/{m + 1), because a pos- 
sible way to use p as a teleportation channel is by using 
a twirled version of an optimal state P{p) and the stan- 
dard teleportation scheme, which produces a depolarizing 
teleportation channel with d = mD{P+, P{p))/ {m + 1) 
(recall eq. (^6[)). On the other hand d„iin{p) is at least 
mA{P+,p)/{m + 1). Indeed, take an optimal teleporta- 
tion scheme employing the state p and LOCC only. It will 
produce some optimal teleportation channel A, that we 
can turn into a depolarizing channel without increasing 
d(Ap'^P) = dminip)- Then we can send half of P+ through 
the channel to obtain a noisy singlet pp that satisfies eq. 
(p6|). The desired inequality follows then from the fact 
that D{P+,pp) > V{P+,Popt{p)). 



Therefore, using ( p9| ) we can announce the follow- 
ing upper bound to the optimal teleportation distance 
dmin{p) achievable with state p and LOCC 

Proposition 6: 

dm^„(p) > (to- l + 2AA(p)). (38) 

The two results of this section can also be derived in 
terms of fidelities (the so-called singlet and channel fideli- 
ties, see for instance ^^). The upper bounds one obtains 
read, 

Fopt = max(<I>+|P(p)|$+) < 1±^:^; (39) 
P m 

foptifi) ^ max / dmmmm < ^'^(•^^^) + ^) . (4o) 

Ap J m + 1 

IV. UPPER BOUND TO DISTILLATION RATES 

We now move to consider a second application of the 
previous measures, namely a bound on the asymptotic 
distillability of a mixed state p in terms of E_\f{p). 

The distillation rate of a bipartite state p is the best 
rate at which we can extract near perfect singlet states 
from multiple copies of the state by means of LOCC. The 
asymptotic (in the number of copies) distillation rate is 
the so-called entanglement of distillation Ed{p) one 
of the fundamental measures of entanglement. In this 
section we will show that the logarithmic negativity £V 
is always at least as great as the entanglement of distil- 
lation E'j~i{p), where e denotes the degree of imperfection 
allowed in the distilled singlets. 

Let T denote a maximally entangled state of two 
qubits, and consider, for some number of copies of 
p, the best approximation to copies of T that can be 
obtained from by means of LOCC. Similarly to the 
previous section we define 

A(T®"°,p^"°) =inf ||T®"° -P(p®"°)||i , (41) 

where P runs over all deterministic protocols built from 
LOCC. We say that c is an achievable distillation rate 
for p, if for any sequences na,ma — cxd of integers such 
that limsup^ (nc/ma) < c we have 

lim A(T^"°,p®"°) = . (42) 

a 

The distillable entanglement Eo (p) corresponds then to 
the supremum of all achievable distillation rates. Several 
variants of this definition are available in the literature, 
which are however, equivalent to the one given here. In 
particular, we may replace 'A ^ 0' by a 'fidelity^ 1', 
and we may consider selective protocols, in which oper- 
ations produce variable numbers of output systems on 
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the same input, and the expected rate is optimized. Of 
course, restricting the amount of classical communication 
between Alice and Bob will in in general change the rate. 

The above definition requires that the errors go to zero, 
but in many applications one can live with a small but 
finite error level. Therefore we introduce Ejj{p), the dis- 
tillable entanglement at error level e, which is defined 
exactly as above, but is replaced by 

lim sup A(T®"° , ) < e . (43) 

a 

Of course, E^{p) = Ed{p), and e ^ E'^oip) ^ i^O'^" 
decreasing function. The main result of this section is 
the following bound. 
Proposition 7: 

E},{p) < E^{p) , (44) 

for aU < e < 1. 

Proof: The only property of LOCC operations used in 
the proof is that for any such operation P, there is an- 
other, P' such that P(p)^-* = P'ip'^^)- We denote by 
the maximally entangled state on a pair of d-dimensional 
spaces. Then, as shown below, we have ||TJ"*||i = d. In 
some sense this is the worst case: for general hermitian 
operators we have < c?||yl|ji. 

Now suppose that P is the transformation for which 
the infimum ( pl| ) for A{T d,p) is attained. Then 

||/"||i>||P'(/-^)||i = ||P(p)'^^||i, (45) 

where the first estimate holds, because P', as a bona fide 
LOCC operation, does not increase the trace norm (recall 
the monotonicity of A/'(p)). On the other hand, 

||P(p)^^||i>||Tj-||i-||(T,-P(p))^-|li 

>d-dA{Td,p). (46) 

Taking the logarithm, we find 

EMp) > log2(d) +log2(l - A(Td,p)) . (47) 

Now let Ua , TOq, be diverging integer sequences as in the 
definition of achievable rate c. Then, using the additivity 
of E^, and the last inequality with d — 2"", we find 

EMp) = — i?A^(p®™°) 
rria 

> -^(n,+log2(l-A(T®"°,p®™°))) . 

We now go to the limit superior with respect to a, ob- 
serving that the error A is uniformly bounded away from 
1, and TUa — > oo. Hence E^{p) > c for every achievable 
rate c, which concludes the proof. 



V. EXPLICIT EXAMPLES 

In this section we display explicit expressions for the 
negativity for some particular classes of bipartite states, 
namely for arbitrary pure states, for mixed states with 
a high degree of symmetry, and finally also for Gaussian 
states of a light field. 

A. Pure states 

All entanglement measures based on asymptotic distil- 
lation and dilution of pure-state entanglement, in partic- 
ular the entanglement of formation Ep and the distillable 
entanglement Ejj but also the relative entropy of en- 
tanglement [|| agree on pure states, where they give the 
von Neumann entropy of the restricted states. Negativ- 
ity gives a larger value: Let p = | | be a pure state, 
and write the wave vector in its Schmidt decomposition 
<i> — J2a ^a^'a ^ i where Cq > are the Schmidt coeffi- 
cients of $, and the ei''' are suitable orthonormal basis. 
Then we get the following result. 

Proposition 8: 

a 

This is precisely Afss/'^, i-e. half of the robustness of 
entanglement, as computed in [ p^ . 

Proof: Introducing the operators "flip" Fe^ (8) = (g) 
e^, and C = J2a'^a\e'a){ea\, and a similar C" for the 
second tensor factor, we find 

= F (C ® C") . (49) 

From the trace norm ||A||i = tiVX^X we may omit uni- 
tary factors like F, so the trace norm is equal to the trace 
of the positive operator (C (g) C"), namely {J2^ Ca) . 

Since £V is an upper bound on the distillation rate, 
and that rate is known to be E{p), the von Neumann 
entropy of the restricted state, we know that Ej\f{p) > 
E(p). But of course, we can get this more directly: using 
the concavity of the logarithm, we get 

Eip) = 2 J^cllog^i-) <2\ogJj2ca) ^Eu{p) . 

a. a 

(50) 

This derivation also allows the characterization of the 
cases of equality: Since the logarithm is strictly concave, 
we get equality if and only if all non-zero are equal. 
Hence equality for pure states holds exactly for maxi- 
mally entangled states (which may have been expanded 
by zeros to live on a larger Hilbert space). 
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B. States with symmetry 



C. Gaussian states 



All entanglement measures can be computed more eas- 
ily for states that are invariant under some large group 
of local unitary transformations [p2|j29| ] . The negativity 
is no exception. The main gain from local symmetries 
is that the partial transpose lies in a low dimensional 
algebra, and is hence easily diagonalized. For this back- 
ground we refer to . But often a direct computation 
is just as easy. 

Consider, for example, the states p on C*, which 
commute with all unitaries of the form U ®U , where U 
is real orthogonal. These can be written as 

p = ad|$+)($+|+5F + cI, (51) 

where |$+) = (X]a=i 1"^ ® Q^))/"^ is again the stan- 
dard maximally entangled vector, and a, 6, c are suitable 
real coefficients. This family includes both the so-called 
Werner states |2l] with a — and, with 5 = 0, the 
so-called isotropic states Q (or noisy singlets, compare 
Eq. ( p7| ) above). The three operators in this expansion 
commute, so all operators of the form (^ can be diag- 
onalized simultaneously, with spectral projections 

Po = |$+)($+| 

Pi = (I-F)/2 

P2 = (I + F)/2- |$+)($+| . 

We parameterize the states of the form (|5^) by the two 
expectation values / — (itr(/9|<i>+)($+|) and g — tr(/9F), 
the third parameter for determining a, 6, c being given 
by the normalization. Then the states correspond to the 
triangle G < f < d, -I < g < 1, f < d{l + g) /2. 

Since partial transposition simply swaps the operators 
F and (i|<I>+)($+|, leaving I unchanged, we can apply the 
same method to compute the trace norm of the partial 
transpose, and hence N{p). Explicitly, we get 

AA(p) = \\l-f\ + \\l + f- 2g/d\ + ^\g/d\ - 1 (52) 

It turns out that in this class of states the Peres- 
Horodecki separability criterion holds (in spite of arbi- 
trary dimension d), i.e., the set of ppt-states is the same 
as the set of separable states, and in the parameteri- 
zation chosen equal to the square f,g E [0,1]. Hence 
Afssip) — J^ppts{p)- Evaluating a simple variational ex- 
pression we get 

A/;4p) = ^max{|2/-l|-l,|2.g-l|-l,0} . (53) 



Gaussian states frequently occur in applications in 
quantum optics, where they describe the light field. Al- 
ice's and Bob's systems are then described by a certain 
number of canonical degrees of freedom such as field 
quadratures of suitable modes. However, the same for- 
malism applies when the canonical operators are posi- 
tions and momenta of a certain number of harmonic os- 
cillators. Gaussian states are then defined as those (pos- 
sibly mixed) states with Gaussian Wigner function. 

For simplicity we denote the full collection of canoni- 
cal operators by a = 1, . . . , 2n, where Alice holds ua 
oscillators and Bob holds n b , and n = ua + tlb ■ These 
are either position or momentum operators, whose com- 
mutation relations are of the form 

[Rc.,Rp]=ia^l,t, (54) 

with an antisymmetric scalar matrix a, called the sym- 
plectic matrix, which has the block matrix decomposition 














(55) 



with respect to a decomposition of the set of indices into 
Alice's and Bob's. This form expresses the fact that all 
variables of Alice commute with all of Bob's. 

A Gaussian state is determined by its first two mo- 
ments ttIq, — iv{pRa) and 

1 

lap = tr{pRaRf3) ~ -(Taf3 , (56) 

where the subtraction is chosen as the antisymmetric 
part of tT^pRaRjs), which is fixed by the commutation 
relations, independently of the state. 7 is then a real 
symmetric matrix. Since the mean ma can be made zero 
by a local unitary transformation (a translation in phase 
space) it is irrelevant for entanglement, and we choose it 
to be zero. The second moment 7 is then the same as the 
covariance matrix of the state. The uncertainty relation, 
the universal lower bound on variances, is then expressed 
as the positive definiteness of ja/i + (*/2)o'q-/3. It will be 
crucial for later that for a classical Gaussian it is only 
necessary that 7 itself is positive definite. Hence we can 
have non-positive operators, whose Wigner function is an 
ordinary, if somewhat sharply peaked Gaussian. 

In order to compute the trace norm of such an op- 
erator or, more generally, to compute the spectrum or 
other characteristics not depending on the Alice — Bob 
partition of the system, we can bring 7 into a standard 
form by a process known as symplectic diagonalization 
or normal mode decomposition. This means choosing a 
suitable canonical linear transformation (i.e., a transfor- 
mation leaving the symplectic form a invariant), which 
can be implemented on the Hilbert space level by unitary 
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operators (known as the metaplectic representation) . As- 
suming a to be in standard form, i.e., block diagonal 

with n 2 X 2 blocks of the form ^ g ) this results in 

a diagonal 7, with equal eigenvalues for each block, i.e. 
7 = diag(ci,ci,C2,C2, . . . ,c„,c„). We cah (ci,C2, . . . ,c„) 
the symplectic spectrum of 7. The fast way to com- 
pute it is via the eigenvalues of the matrix (T~^7, which 
are ±jci, . . . ,±ic„. At the Hilbert space level the nor- 
mal mode decomposition transforms the state into a ten- 
sor product of independent harmonic oscillators, each of 
which is in a thermal oscillator state, the temperature 
being a function of Cq. The smallest value allowed by un- 
certainty is Cq, = 1/2, which gives the oscillator ground 
state. 

In this context transposition is best identified with time 
reversal. Indeed, in the general scheme we can choose a 
basis in which the transpose is computed, but all these 
choices are equivalent via a local unitary transforma- 
tion. In this case we choose the position representation, 
in which transposition is the same as reversing all mo- 
menta, keeping all positions, and to lift this to prod- 
ucts by observing that transposition reverses operator 
products. Partial transposition Ta is completely analo- 
gous. Only in this case just Alice's momenta are reversed 
and Bob's arc left unchanged. On the level of Wigner 
functions and covariance matrices of Gaussians, we just 
have to apply the corresponding linear transformations 
on phase space. That is, 7"^^ , the covariance matrix of 
the partial transpose of a Gaussian state with covariance 
matrix 7 is constructed by multiplying by —1 all matrix 
elements, which connect one of Alice's momenta to either 
a position, or a degree of freedom belonging to Bob, and 
leaving all other matrix elements unchanged. 

The point is, of course, that while this transformation 
preserves the positive definiteness of 7, and hence we get 
another Gaussian Wigner function, it does not respect 
the uncertainty relation, so the partially transposed op- 
erator may fail to be positive. However, the whole formal- 
ism of normal mode decomposition for 7 works exactly as 
before: we get a representation of the partial transpose 
as a tensor product of (not necessarily positive) trace 
class operators. The trace norm of this operator is just 
the product of the trace norms, so we have completely 
reduced the computation of the trace norm to the single 
mode case. To summarize the results so far: 

Let p be a Gaussian density operator with covariance 
matrix 7, and let (ci, . . . ,c„) be the symplectic spectrum 
of 7^^ . Then 



(57) 



where F{c) — log2||pc||i, md pc is the operator whose 
Wigner function is a Gaussian with covariance diag(c, c). 
Of course, the function F vanishes for c > (1/2). It is 



easily determined by looking at Gaussian states for oscil- 
lators as the temperature states of the oscillator. With 
z = e~^, and \n) the nth eigenstate of the oscillator, a 
general Gaussian is of the form 



p^{l-z)J2z- \n){n\ 

n=0 



(58) 



where z > corresponds to density operators, and — 1 < 
z < to Gaussians whose Wigner functions have sub- 
Heisenberg variance. Then we get 

-1 



IIHI 



(l-z)(l 
tripP^) = 



(59) 



tr(pl(p2 



+ 7;) 



Solving Eq. 

F(c) = 



|) for z and substituting into (| 

for 2c > 1 
-log2(2c) for 2c < 1 



(60) 
t) we find 
(61) 



Together with ( |57| ) and the process of normal mode de- 
composition this is an efficient procedure for determining 

In the simplest case of one oscillator each for Alice and 
Bob we may go even further, by expressing the symplectic 
spectrum of the partial transpose directly in terms of the 
covariance matrix. Suppose that 



7 



A 






B ^ 



(62) 



with 2x2 matrices A,B,C. Then, as shown in [ pi| , 
the numbers det^, detB, detC , and det7 are a com- 
plete set of invariants for 7 with respect to local symplec- 
tic transformations. Moreover, when passing from 7 to 
7"^-^ only det C changes sign, and the others remain un- 
changed. Pure states are characterized by the conditions 
det7 = 1/16, and det A -I- det B -|- 2 det C = 1/2, and can 
be brought into the normal form 



/a 




c 

Vo 



\ 



(63) 



where = -I- 1/4. 

Coming back to the general case of (^), the charac- 
teristic equation of a^^j"^^ , whose solutions are the ±Ca, 
takes the form 

+ (det A + dot B - 2 det C)^^ + det 7 = (64) 

Together with Eq. (|5^) this amounts to an explicit for- 
mula. For the particular case of a pure state we find 



E^ip) = -21og2(Va-l/2- V« + l/2) 
which is readily seen to agree with Eq. (El). 



(65) 
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VI. MULTIPARTITE SYSTEMS 

As argued in the introduction, a computable measure 
of entanglement for bipartite mixed states is also very 
convenient for the quantification of multipartite entan- 
glement. In this section we describe a whole set of com- 
putable parameters related to the negativity that can be 
associated to a multipartite state to make quantitative 
statements about its entanglement. 



A. Multipartite negativities 

Consider a quantum system consisting of, say, three 
parts, associated to Alice, Bob and Charlie. Let pABC 
be the (either pure or mixed) state of the system. A pos- 
sible way to classify the entanglement properties of such 
a state is by looking at the different bipartite splittings 
1^ of the system. 

First, we can join two of the three parts, say those 
of Alice and Bob, and compute the sum of negative 
eigenvalues of p'abc^ J^{AB)-c[pabc)- This is auto- 
matically an entanglement monotone ||33| , which quanti- 
fies the strength of quantum correlations between Char- 
lie and the other two parties. Similarly, the negativi- 



ties N{^ac)-b{pabc) and N{^bc)-a{pabc) are two other 
monotonic functions under LOCC with analogous mean- 
ing. 

We can also consider the entanglement properties of 
two-party reduced density matrices. Suppose, for in- 
stance, that Charlie decides not to cooperate with the 
two other parties in the manipulation of the tripartite 
system according to LOCC. Alice and Bob's effective 
density matrix, gab =Trc[pAsc], may still retain some 
of the original entanglement. The negativity of ctab, 
J^A-B:g:{PABc)i Can be used to quantify this residual en- 
tanglement. Analogous quantities can be used to quan- 
tify the entanglement of <tac and (Jbc- 

Thus, altogether we have obtained 6 computable func- 
tions to quantify the entanglement of any state of a tri- 
partite system. In a four-partite setting the number of 
possible splittings is much greater (see |32 for a more 
detailed description), and thus, we obtain up to 26 in- 
equivalent measures, namely: {i) Na-bcd{pabcd) and 
the corresponding permutations, i.e. 4 inequivalent mea- 
sures; ill) Nab-cd{pabcd) and permutations (4 mea- 
sures); (iii) Ma-bc;p{pabcd) and permutations (12 
measures); {iv) J^a-B;0p{pabcd) and permutations (6 
measures). 



B. Hierarchy 

Notice that although all these measures are inde- 
pendent functions of the multipartite state, there is a 



strength hierarchy between them when corresponding to 
related bipartite splittings with different number of par- 
ties. In the four-party case we have that, for instance. 



AC 



A-BCD 



A-BC;l 



A--B;^P, 



(66) 



which follows from the fact that to trace out a part of a lo- 
cal system is an operation of the set LOCC, under which 
the negativity can only decrease. Of course, the same 
inequalities hold for the corresponding logarithmic nega- 
tivities, and thus also for the several bounds on distilla- 
bility — of different kinds of multipartite entanglement — 
implied by the later. 

It should be noted, however, that in this way one can 
quantify only some aspects of multi-particle entangle- 
ment: there are tripartite states which are separable with 
respect to every splitting of the system, but are neverthe- 
less not a convex combination of triple tensor products 
of density operators States which have positive 

partial transpose with respect to every subsystem satisfy 
a large class of Bell inequalities . 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we have presented a computable measure 
of entanglement for bipartite mixed states, the negativ- 
ity Af{p), which we have proved not to increase under 
LOCC. Although it lacks a direct physical interpreta- 
tion, we have shown that it bounds two relevant quanti- 
ties characterizing the entanglement of mixed states: the 
channel capacity and the distillable entanglement Ejy. 

Ideally, quantum correlations would be best quantified 
by measures with a given physical meaning. Which mea- 
sure to be used, would depend on which question we want 
to answer. For instance, if we want to know how much 
pure-state entanglement the parties can extract from (in- 
finitely) many copies of the state p, then the proper mea- 
sure to be used is the entanglement of distillation Eu{p)- 

In practice, however, the value of these measures is 
not known. Recent studies of entangled systems, such as 
those of entangled chains, entanglement molecules, entan- 
gled rings, entangled Heisenherg models and cluster states 
p6| , which are A^-qubit systems in some global entangled 
state, are focused on the two-qubit quantum correlations 
associated to the global state, as measured by the en- 
tanglement of formation (or the related concurrence) [|| . 
This choice of measure of entanglement is somehow arbi- 
trary — where not simply forced by the lack of an alter- 
native measure that can also be computed for two-qubit 
mixed states — , because it does not reflect anyway the 
entanglement cost of formation of the A^-qubit state. We 
envisage that in these and similar contexts it will pay off 
to use a computable entanglement measure like the neg- 
ativity, whose evaluation is not restricted to two-qubit 
mixed states. The negativity will allow, for instance, to 
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generalize the previous investigations to analogous con- 
structions with Mevel systems {I > 2) instead of qubits, 
also to analyze quantitatively the entanglement between 
subsets of these /-level systems. 

Finally, in a similar way as the negativity has recently 
played a role in proving the irreversibility of asymptotic 
local manipulation of bipartite mixed-state entanglement 
| p6[ , we hope that this computable measure will also be 
a useful tool to answer other fundamental questions of 
entanglement theory. 
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